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GENERAL THEORY OF LINEAR DIFFERENCE EQUATIONS* 

BY 

GEORGE D. BIEKHOFF 

Introduction. 

The theory of linear difference equations with rational coefficients was in a 
very backward state until Poincar6 f in 1882 developed the notion of asymp- 
totic representation, and its application to this branch of mathematics. Further 
important progress in this direction has been made only very recently, notably 
by Galbrun,J who by means of the Laplace transformation has investigated 
the properties of a set of fundamental solutions in the entire plane of the com- 
plex variable x . 

The aim of the present paper is first to study the nature of these solutions to 
which I have been led by means of direct methods, and secondly to show that 
there exists a purely Eiemannian theory of linear difference equations. In par- 
ticular certain rational functions of e^"^""' are shown to play a part like that 
of the monodromic group constants of an ordinary linear differential equation. 

To the best of my knowledge, the importance of the functional standpoint in 
the field of difference equations was emphasized first by Van Vleck in an 
inspiring series of lectures given at the University of Wisconsin in the spring 
of 1909, in which he conjectured the existence of sets of solutions analytic on 
either the left or the right side of the complex plane. 

On account of the extreme simplicity of the matrix notation, I have found it 
convenient to deal with a linear difference system of n equations of the first 
order 

(1) 9^{^+^)=i%i^)9j{^^)^ (i = l,2,...,n), 

rather than with a single equation of the nth order. In these equations the 
functions a..(x) are taken to be rational functions of x with a pole at a; = oo of 
order /x, at most, so that one has 

(2) a,.{x) = a^.x'^ + ^'jcC^-' + • • •, {\x\>B). 

* Presented to the Society, September 6, 1910. 

t American Journal of Mathematics, vol. 7 (1885), pp. 203-258. 
JComptes Rendus, vol. 148 (1909), pp. 905-907. 
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The sole other restriction to be imposed on these functions below will, in par- 
ticular, be satisfied if the determinantal characteristic equation 

(3) |«,-8,./>J = 

[where S. = 0, i '^j; B^. = 1] has n distinct roots which are not zero. If the 
n sets of functions 

are n linearly independent solutions of (1), the array of elements g^{x) forms a 
matrix solution G{x). The array of elements a.. (a; ) forms a second matrix 
A[x), and the n^ equations which the n solutions satisfy can be combined into 
a single matrix equation 

(4) G{x + l)=A{x)G{x). 

If X be replaced by a; — 1 , this may be written 

(4)' G(x-l) = A-\x-l)G{x). 

The most general matrix solution of (4) is furnished by 

(5) H{x) = G(x)P(x), 

in which G(x) is any particular matrix solution, and P{x) is any arbitrary 
matrix of periodic functions of period 1 whose determinant does not vanish 
identically. 

The equation (4) admits of two symbolic solutions 

G{x) = A{x-l)A{x-2)..., 

G(x)= A-\x)A-\x + I) . . ., 

which converge to limit matrices only in particular cases. This convergence is 
secured to a sufficient extent by a suitable modification, made in § 1, as follows : 
The system (1) is known to be satisfied formally by taking for ^j( a;), ^2(0!), 
■ • •,g^[x) certain formal series, one set of them for each of the n roots p of the 
characteristic equation (3). The matrix formed from these sets taken in a cer- 
tain order is denoted by S{x), and the modified matrix obtained by breaking 
off these elements at the ^th term is called T(^x). The sequences of matrices 

A(x-l)A{x-2).-.A(x-m)T(x-m), 

(m=l,2, •••), 
A-\x)A-\x + l)-'-A-\x + m-l)T{x + m), 

form the desired modification of the two symbolic solutions. It is demonstrated 
that all the determinants formed from the first \ columns and any X rows of the 
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first of these matrices, ( X = 1 , 2 , • • • , n ) , or from the last X columns and any 
X, rows of the second of these matrices, converge to limit functions independent 
of A; as m becomes infinite. By the aid of these functions and by a process of 
summation based on a contour integral suggested by the one which Carmichael 
has used, two remarkable particular matrix solutions G(x) and H{x) are 
obtained in §§ 2, 3, 4. These are termed the first and second principal matrix 
solutions respectively. The elements of G(x) are analytic, save for poles at 
points a -\-l, (7 + 2, •••, where o- is a pole of one of the elements of A(x), 
and have the asymptotic form of the corresponding elements of jS{x) in any 
left half plane. The elements of Hi^x) are analytic save for poles at points 
<T — l,cr — 2, ••., where cr is a pole of one of the elements of A~^ {x — 1 ) , and 
have the asymptotic form of the corresponding elements of /^(a;) in any right 
half plane. The principal matrix solutions G(x) and Il{ x ) are uniquely deter- 
mined by this asymptotic property. A relation (5) will also, of course, subsist 
between G{x) and II{x'). 

Conversely if G{x) and II{x) are asymptotically represented by the elements 
of a matrix 8{x) even to terms of the first order and are related by means of 
a matrix P{x) of periodic functions, as in (5), these matrices G{x) and Il(x) 
will form the first and second principal matrix solutions of a system (4) in which 
the elements of A(x) are rational (see § 7). 

It is proved in § 5 that the elements of P(x) are rational functions of e'^'^^'^" . 
The structure of these functions is determined with particular reference to the 
case when the functions aJx) are polynomials of degree /i, to which case the 
more general one is reduced. 

These results lead in § 6 to the determination of the nature of the elements 
of G(x) and If( x) in the remaining right and left half plane respectively ; the 
asymptotic form is found to alter along certain critical rays. A complete 
description of the principal matrix solutions is thus obtained. 

In § 7 it is shown that the number of characteristic constants involved in the 
complete characterization of G(x) and H{^x) and their connection by means 
of i-'(x) is precisely equal to the number of arbitrary constants in the coefiicients 
a^.{^x) assumed to be polynomials, and that the characteristic constants are inde- 
pendent. Thus I have ventured to propose the following as a fundamental 
problem of linear difference equations : to determine whether there exists a linear 
difference system with prescribed characteristic constants. 

It would be interesting to know the exact relation between the characteristic 
constants and the monodromic group constants of the differential equation yielded 
by the Laplace transformation. 

The method contained in the present paper is of wide generality and applies 
if the coefficients a^.(^x) have the character of rational functions only at a; = oo, 

* These Transactions, vol. 12 (1911), pp. 99-134. 
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in all cases where there exist n formal solutions. Carmichael (loe. eit.) has 
studied this ease to a considerable extent ; in particular he has proved the exist- 
ence of analytic solutions and has partly discussed their asymptotic form in the 
direction of the axis of reals. 

§ 1. The Determinant Limits. 

It is apparent that the given linear system (1), when written in the matrix 
notation (4), (4)', possesses the two symbolic solutions 

A{x-l)A(x-2) ..., A-\x)A-\x-{-l)---. 

Ihese will converge to limit matrices in certain special cases, for example if 
ff..(cc) has the form 

in this case the above symbolic solutions will yield the complete solution of (1). 
In general however these expressions must be properly modified in order to be 
made to lead indirectly to a complete solution of (1). 

Let us first make precise the limitation which is to be made on the coefficients 
«..(»;) in equations (1). If the roots of the characteristic equation (3) are dis- 
tinct, and if none of them are zero, there will exist precisely n sets of series 



s,.{x) = x-'^ip.e-'^ fx-' \ Sy + ^ + 



s' 



(6) «2,-(»') = ^'^'(Pye~'')^»'''M«2y + 5'+ •••}' 0- = l,2,.-.,»), 



«(i) 



s^.{x) = x^^ipjc-- yx'-^ Y„ + ^ + 

each of which constitutes a formal solution of (1). The proof of this is a matter 
of direct reckoning, and the fact will be taken for granted.* The constants 
/3j, p^, • • • , p^ are the roots of the characteristic equation, and furthermore we 
have the determinant 

(T) c?=|s,|4=0. 

* The formal solutions of a single linear difference equation of the nth order have been fre- 
quently used (for example see Horn, Crelle's Journal, vol. 138 (1910), pp. 159-191) and 
admit of easy extension to a linear system. The formal solutions may be obtained by direct 
substitution if the roots of the characteristic equation are distinct and different from zero ; the 
series expansion 






2x 

enables one to make a comparison of coeflBcients and thus to determine the constants of the 
formal solutions. 



a!" 



.n»a; / 
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The restrictive hypothesis which we shall make is not that Pi, p^i • • •■> Pn '^'''^ 
distinct hut only that there exist n formal solutions (6) with the speckled prop- 
erty (7). The quantities p^, p^^ — ? />„ must of course still be the roots of the 
characteristic equation, and none of them can be zero. We shall choose our sub- 
scripts so that 

(8) \pA^\PA = ---=\Pn\- 

Without loss of real generality we may add the further condition that the 
linear homogeneous difference equation of the wth order at most which ^,( a;) 
satisfies, and which is obtained from (1) by elimination, is actually of the nth 
order. 

The elements s^^x) form a matrix S{x) which is a formal matrix solution of 
the equation (4) 

S(x + l) = A{x)S(x). 

Let us compute the elements s'.j(x) of the inverse matrix S~^(x). The deter- 
minant I yS'(x)| may be written 

"ipip2 • • • p^e-'"^rx'-^+'-^+ ■■■+'■•' I <z -^ — -K . . . I . 

The element in the *th row and Jth column of the inverse matrix is the quotient 
of the cofactor of the element in the jth row and the itla column of this deter- 
minant by the determinant itself, and is therefore given by a formal series 

(9) s:^{x) = x-'^^p.e-'^ y^x-' I ..^^ + ^ + . . . [ . 

Let T[x) denote any matrix obtained by breaking off the elements of S(x) 
at the kth term, or more generally by replacing s^(x) by convergent series t.j(x) 
of the same form as s^.ix) and having the same first k terms. Write then 

(10) T(x+l) = B{x)T{x). 

The matrix B(x), as thus defined, is clearly a matrix of functions J>„(^x) each 
of which has an expansion like a.. (a;) at x = oo with precisely the same first k 
terms. For if we compare the formulas 

B{x)=T(x + l)T-'(x) and A (x) = S {x + 1) S-\x) , 

we observe first that the element t'^.^x) of T~^(^x) is given by (9) to the first 
k terms, and that in consequence the right hand members of both the equations 
have the same first k terms. In consequence we may put 

(11) A{x) = B{x)ll+x->'C(x)}, 

1 being the identity matrix, and the elements of C( a;) being analytic at «;= oo . 
It is to be noted that the determinant formed from the principal coefficients of 

Traas. Am. Math. Soc. 17 
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A(x) and .B(x), namely | a;.|, cannot be zero, for that implies that zero is a 
root of the characteristic equation. 

In this paper it is necessary to distinguish sharply between two kinds of 
asymptotic representation of a function by a series of the form (6), say 

s(x) = x>'^{pe-'^fx''l s + —-+•• • ). 

Let g(x) be the given function ; if for each Jc the difference 

g{x)x->^-{pe->^)-''x-- - (^s + ^ + . . . + ^ j 

becomes uniformly small of order x~* for lim x = oo in a certain region, we say 
that g (x) is asymptotically represented by s(x) in that region, with respect to 
X ; if on the other hand the difference is uniformly small of order v~'' (where 
x — u+ l/ — Iv, u, V being real quantities) for lim tj = oo in the region, we 
say that the representation is with respect to v. In this connection two 
important inequalities are 



(12) 



(^x-vf\^2 



for X in the left half of the complex plane, and 
(13) -^ ^ 



(x — vy 



2 rj^ tt] 



(^S2), 



(^S2), 



for X in the right half of the complex plane. To prove the first of these we note 
that when x lies in the left half plane u is negative, so that 

\x — v\'' = (u — vy + V' > u^ + v^ + 1^ , 

whence the sum of the given series is less than 

dv TT 



» 1 T" 1 



{\xf+j^)~2\x\''-'' 
Again, for x in the right half plane, our series may be written 

1 +^1 1+^1 



-( [(^u-vf + v^y 



fi< E 



[(u-vf + v^y 



1 +00 

ii, ^ I ., I*— 2 2-1 



■^^{u-vf + v'' 



If each quantity | m — z^ | be replaced by the last of the integers 0,1,2, ... 
which does not exceed it, the terms in the last series are not decreased, and hence 
one has 



1 



+» 

E 



{u 



1 ^ 2 " 1 2 r 1 r" dv 
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But the last integral reduces to 7r/2 | v |, and thus the second inequality holds 
also. 

We may now proceed to the existence theorem : 

Theorem I. Form the two sequences of matrices 

P^{x) = A{x-\)A{x-2) . . . A{x- m)T{x - m), 
(14) (m=l,2, •••). 

Q^{x) = A-'-{x)A-\x + l)-.-A-\x + m-l)T{x+m\ 

Every determinant formed from, the first X columns (\= 1 , 2, • • ., re) and the 
ith, jth, • • • , Ith rows of P^ {^)\.Qmi^)'\ converges, for h sufficiently large, to a 
definite limit function w^-...i(a') ['W^... «(«')] , independent of h, as m becomes 
infinite. This function is analytic throughout the entire x-plane excepting at 
points any number of units to the right (left) of poles of the elements of 
A(^x) [^~'(x — 1)] and at these points may have a pole. The asymptotic 
form of '>hj_i{x)\^Vij___i{x)'\ in the vicinity of x = oo is given by the corre- 
sponding determinant %...i(x) of S{x), with respect to x in any left {right) 
half of the plane, and with respect to v in any right (left) half of the plane. 

Proof. We shall first prove the theorem in so far as it relates to the case 
\ = 1 , and then show how the cases X = 2 , 3 , • • •, n may be reduced to this 
case. 

It is necessary first to show that the elements of the first column of P^(x) 
converge to definite limit functions u^{x), u^{x), • • ., w^(a;) as m becomes 
infinite. 

The matrix P^{x) may be written 

P^{x)=T{x)P^{x) 

= T{x){lT-Xx)A{x-l)T{x-l)] [T-\x-l)A(x-2)T{x-2)'] 

... [T-\x-m+ l)A(x — m.)T{x-m)]}. 

It is sufficient to prove that the elements of the first column of P^{x) will con- 
verge, inasmuch as the elements of the first column of P^{x) arise from com- 
bination of the rows of T{x) with the first column of P^{x). The matrix 
P^{x) is a product of matrices, the type of each of which is 

T-\x + l)A{x)T{x). 

By (10) and (1 1) this reduces to 

l + ~T~\x)C{x)T{x). 

Let us examine the second term. The element in the ith row and jih column 
of T~^{x) will be given by an expression like that for s'^Ax) in (9) with the 
same first h terms; the elements of C{x') are analytic at a; = oo , the element 
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in the ztli row and^'th column of T(x) will be given by an expression like that 
for s^,(x) and with the same first k terms. The element in the ^th row and^'th 
column of the product T~^{x) C{x) T(x) has therefore the form 



{^J"'- 



times a function analytic at x = oo . 



Choose first k so large that it exceeds every difference r. — r^ in absolute value 
by more than an assigned integer cZ = 2 , and we see that the second term may 
be written 

(16) ^-,@{x) with elements 0..{x) ^ (^^^^ cf,..{x). 

The function 4'f{x) is finite in the vicinity of a; = oo, at least if x be restricted 
not to pass the positive axis of reals. 

By means of this notation, P^^ ( a; ) reduces to the form 

or more fully. 



(16) 



h<h 



+ z 



ttA^-hYi' 



^^e{x-k^)€>{x-k.^)+ .. 



It is the elements of the first column of this matrix which nmst be shown to 
converge. According to the notation of (15), the ith element is 



(17) 






But (f>i.{x) is limited in the vicinity of :<; = oo and is less than some constant 
JbTfor I a; I > -ffj . Therefore if all the points a; — 1 , a; — 2 , • • • are exterior to 
the circle \x\ ^ li^, which is true if x lies in D of fig. 1, the typical element in 



\x[^Bi 



the(^+ l)th termof (17), 



Fig. 1. 
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{x — k^f {x — h^f {x 

X 4>,,{x-k^)4>^^{x-h^)...4>^^{x-k;), 

is less in absolute value than the quantity 

(Py\ M^ 

\pj {x-k,y...{x-k,y 

This inequality depends on the fact that we have 



fP_r_Y'"'(P_.Y'\..(P-lY~'" 



m(ssi^r-(s^"'>m 



as a consequence of (8), and the fact that k^, k^, • • • , k^^ are an increasing set of 
positive integers. Hence each term of (17) does not exceed the corresponding 
term of 



(18) S,+ 



m 



M 



■ k 



T4 + 



m 






nM'- 



i> K* rvt , A" 



r^ + 



in absolute value. 

Allow m to approach infinity. The expression (17) becomes a multiple series. 
Since the sum of the terms after the first in (18) approaches the limit 



(19) 



a;)1l[(>-i^)(-^)-]-l' 



it is obvious that the elements of the first column of P^^^(x) and therefore of 
J-'^^(x) converge absolutely and uniformly to an analytic function in the vicinity 
of points in D. Moreover as we may write 

I>^X^) = A{x-l)A{x-2)...A{x-r)P^^^_^,{x-r), 

and as we may take r so large that x — r — l,x — r— 2, • • ■ lie without the 
circle \x\ = Ji^, the elements of the first column of I^^^(x) will converge uni- 
formly in the vicinity of any point of the plane excepting those only at which 
one of the elements of A(^x — 1), A (x — 2), • • • has a pole. Hence the ith 
element of the first column of P^^(x) converges to a function u.(^x) which is 
analytic except at the points congruent to (any number of units to right or left of) 
the poles of the elements of A(x) on the right. At these points u.(x} may 
have a pole. 

It may easily be made evident that the definition of i<j ( x ) , m^ ( a; ) , . . . , m^ ( x ) 
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is independent of the choice of h and T. For if Jc is any second choice of k 
and T\x) is the corresponding function 7'( x ) , we may nevertheless write 

P:,X'x)=^T{x){[T-\x)A{x-\)T{x-\)^ 

X {T-\x-l)A{x- i)T{x-i)'] ■ . ■ [T-\x-m+l)A{x-m)T\x-m)]. 

If the elements of the first column be expanded in a sum like (17) and m be 
allowed to become infinite, the resulting multiple series for ^l■^{x) will be term 
for term identical with that obtained before. 

It remains to discuss the asymptotic nature of 'u^(x), • • • , u^^(x) for large 
values of |a;| . The limit w.. (x) of the ith element of the first column of P^[x) 
differs from S^^ by a quantity less than (19) in absolute value if x lies in D. 
This quantity is itself less than 

in absolute value. If now any line be drawn parallel to the axis of imaginaries 
and if x be restricted to lie in D to the left of the line, by inequality (12) the 
last factor will be of the form M(x)/x'^~^, where M[x) is limited. Thus we 
may write 

for X in D and to the left of this line. Likewise if a; is in Z) to the right of 
such a line, by inequality (13) we may write 

-.(«') = ^n + (^y'^7S-\ \M\{x)\<M. 
But the relation P^{x) = T{x)P^{x) gives at once 

n 

Ui{x) = J^t.^{x)u^(x) . 

T = l 

If in this equation the first of the above expressions for u^(x),u^(x), • ■ •,u^(x) 
be substituted, we find at once 



^.(«')=<a(»')+|:«..(«')(^y^S- 



The functions ^^.(cc) are the same in form as s..{x) up to the ^th term. Substi- 
tuting in their expansions we obtain 



u,ix) = x-K-^' j(^a + f + • • • j + Z 7^.^[^. + if + • • • jj 
The second term in brackets can clearly be made infinitesimal in x of as high 
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order as may be desired by taking h, and therefore d, large enough. Consequently 
u^(^x) is represented asymptotically by s^^{x) with respect to x to the left of 
any such line and likewise with respect to v to the right of such a line. The 
part of the theorem that relates to the case X = 1 is proved. 

It is clear that if the subscripts 1,2, ■ • •, n had not been so chosen that (8) 
was true, the only modification would be that the elements of any ith column of 
P„{x) for which the quantity p. had a maximum absolute value would 
converge. 

The case X = 2 will next be discussed. This reduces at once to the case 
X= 1 inasmuch as the two-rowed determinants u'-^^^[x) formed J^rom the ele- 
ments of the first two columns of JF*^'"^ (x) which are in the ath and l3th rows 
(a <^) are related to a certain difference system precisely as the elements of 
the first colum7i are related to the original system (1). 

In fact, form the system with general solution 



(20) 



9ij:{^) = 



g'{x) g'Ux) 



Here g[{x), g'^{x), ..., g[Xx) &nAg'[{x), gl{x), 
solutions of (1). One finds by the use of (1) 



{!,i = l, 2, •••, n; i<j). 



g"Xx) are any pair of 



(21) 
where 



<7..,(x + 1) = 23«<,:„,(a')^,,;(a3), 



Accordingly we have ^n(n — 1 ) linear difference equations of precisely the 
type (1) to determine the functions g...[x). It is clear that the sets of formal 
expansions 



ij '■ Jcl 



(x) = 



= x'''%p^p^yx'--+" { ( s.,^ s., - s., s.^) + terms in a;-' , x- 



}. 



h and I fixed, will satisfy (21). There are ^n[n — 1) sets of these, and further- 
more the determinant of the constant terms 



o ., o ., ^^ 6 •, 6 ■, 



(i<j; k<l) 



is precisely equal to d''~^ and is therefore not zero by (7).* Hence the system 
(21) satisfies the restrictive hypothesis made at the beginning of § 1. 

We are now in a position to apply the result already obtained. In matrix 
notation, equations (21) may be expressed as 

G,{x + \) = A,{x)G,{x). 



* Pascal : Die Determinanten, p. 87. 
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Also let T^(x) be the matrix which the formal matrix solution S^i^x) becomes 
when the formal series are replaced by convergent series which agree with them 
to the ^th term. The elements in that column of 

for which the quantity Pj^p^ in the formal solution which appears in the sam,e 
column of S^{x) has a maximum absolute value will converge ; this results 
from an application of the preceding result for X = 1 to the system of Jw ( w — 1 ) 
equations. Equations (8) show that we may take ^' = 1 , ^ = 2 and it is also 
clear that we may take the elements of this column of T^{x) to be 



iij:n\^) ■ 



t,lx) t,,ix) 
t.,{x) t.^{x) 



The elements p':^}n{x) of this column of P^^X^) converge to functions analytic 
everywhere except for poles at points to the right of the poles of the elements of 
A^(^x) and congruent to them. Moreover these functions are represented 
asymptotically by s^j.^^{x) throughout the entire plane, with respect to x to the 
left of any line parallel to the axis of imaginaries, and with respect to v to the 
right of such a line. 

But the sequence pf}i2{x), m = 1, 2, ... is identical with the sequence of 
determinants u''^\x) formed from the elements of the first two columns which 
are in the ith and _/th rows of P^{x). In fact the first member of both 
sequences is <jj:i2(x), and in general pftiii^) is derived from p[p_n{x) just 



as u^+^\x) is derived from u'^'>(x), namely we have 

mm= X (l.rM^ - 1 )pi%.{x - 1 ) , 






H= 



MT+''H. M^+'^H 



'''^+'\x), p^^+'\x) 



Pjl 



E «,,(*- iMT'(x-i), E ^iri^ - i)?gX»' - 1) 



E aj^{x-i)p^:^\x-i), E «,>(«' - i)i'ir(x - 1) 

T = l T=l 

X%:.^(a;-lX'(a'-l)- 



Thus the theorem is true for X = 2 also. Strictly speaking, we have only 
obtained functions M..(a?) for i <.j; but we have obviously such functions for 
all i andj; these may be obtained by means of the equation u^.[x) = — u.^x). 
In like manner the theorem can be demonstrated for the remaining values of X . 
The existence of the functions v^j,,,i{x) and their properties may be estab- 
lished in the same way. The role of right and left in the as-plane will be inter- 
changed and the equation (4)' replaces (4). 
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The functions iiy^ i(x) and i^. ,(a;) will be called the determinant limits 
on account of their origin. 

In the case that the restrictive hypothesis is satisfied and all the roots p of 
the characteristic equation are equal, Theorem I supplies a complete existence 
theorem, every column of elements of P^{x) and Q^{x) converging to an ana- 
lytic solution of (1). These conditions are met if one has 

a..{x) = 8.. + -^+^+ .... 
,j\ J 'J x^ x^ 

In general however Theorem I supplies only one solution directly, and that is 
furnished by 

g^{x)=u^(x),g^{x) = u^{x), ■ . ■, g^{x) = u^{x). 

This set of functions must be a solution since the definition of u^(^x) gives us 

n 

w.(a; + l)=lim X) «ir(*) • • • ^'^^C* — »0^>,i(*~ "0 

m^oo T, - - ., ^^^ 

= ll<-hr{x)u,{x). 

§ 2. The Solutions Associated with the Determinant Limits. 

The functions Mj(x), m^..(x), .•. which have been obtained behave as 
though they were constituent determinants of a matrix of solutions of (2). 

Theorem II. There exist n solutions g^.(^x), ff-2ji^)^ ' ' ' ■> ff,, i^) LKji^)^ 
h^.{x), • • • , /i.^.(cc)](jf = 1 , 2, ■ ■ ■ , n) of (1) over any discrete set of values 
congruent to a given value x^ofx, such that every \-rowed determinant formed 
■from the first \last'\ \-columns of the matrix G {x)\_H {^x)'\ and the ith, jth, 
■ ■ ■,lth rows will equal w_y , ,(cc) [tv.,_,(x)] at these points. The first solu- 
tion ^ii(a?), g^ix), •••, 5'„i(a5)['^i„(*)' ^2«(«')' •••' ^,.„("')] ^« uniquely 
determined, the second solution g^^i^x), g^ii^)-! '''i 9n2i-^)\.^\,n-ii.'-^)i 
^2, M_i ( * ) 5 • • • 5 Ki, n-\ ( * ) ] *'* determined up to an additive term obtained by 
multiplying all the corresponding elements of the first solution by an arbitrary 
function of period 1, a7id so on, any one of these solutions being determined 
up to an additive term linear in the corresponding elements of the preceding 
solutions with arbitrary multipliers of period 1. 

Proof The theorem is true in so far as it relates to X = 1 since Wj(x), 
u^{x), • • • , w^(cc) has been seen to give a solution of (1). We may write then 

g^^{x) = u^{x), gi^{x) = u^{x), ■..,g,^^{x) = u^{x). 

Next we shall demonstrate that there exists some solution g^ii^)^ fftzi^)^ ' ' '■> 
9n2{^) such that for all i andj 
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(22) 



■ u..{x). 



For this purpose we shall recur to the sequence of elements of /^^^(cc). By 
means of the identity 

pf^{x) 2^^{x) p^«>{x) 

p<j.f{x) pf^(x) pi^ix) =0, 

fi;^{^) /4:?(^) p^si^) 

we infer, by developing according to the elements of the third column and letting 
m become infinite, that 

u.j{x)u^{x) — u^,Xx)u.{x) + u.,Xx)u^{x) = (i,j, fc=l, 2, ■ ■ •, «). 

The n quantities g'lsCa'o)' ^22(^^0)' ' ' '■> 9n2(^o) uniquely determine a solution of 
(1) on the set of points congruent to x^. Choose these n quantities so that (22) 
is satisfied for some fixed i such that g.^ ( x^ ) =)= , and for every j . All the 
rest of the equations (22) will then be satisfied, as is evident from the pre- 
ceding equation. This gives us 



ffkii^)- 



9A^) 9j2i^) 
for x = x^, which reduces to 

u.^{x) = 



9n{^) 9i2{^^ 
9ki{^) 9k2{^) 



gji{x) + u.^(x)g,^{x) = 



9jii^) 9j2{^) 

for x = Xf^. Hence g...(x) and u..(x) are equal at x = x^ for all i and^'. 
But they are both constituent elements of a solution of (21), and therefore will 
agree at all congruent points, since only one solution of that equation exists on 
the discrete set of points congruent to x^ whose elements have given values at 
this point. It is necessary to make an exception of points congruent to the 
poles of ^ (a;) or A~^(x — 1), and such points are excluded. 

The method by which this second solution was obtained clearly allows us to 
add to each element g^^i^) a term of the iorm p (x) g ^^(x) , where p{x) is of 
period 1. 

One may now in an entirely analogous way define a solution ^/^(a;), (/'jsC*)' 
• • •, ^^3(0;) on the set of points congruent to x^. The identities to be employed 
in this case are 

«,7A(a')«i(a') - '«i,-;(«')'-^/<(»=) + m„,;(x)m^.(x) -M^.j,(a?)M;(a;) = 0. 

The function g.^ ( x ) is only determined up to an additive term 

p{x)g,^{x) + q(x)g,^{x). 
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By proceeding thus the theorem may be proved in a succession of steps in so far 
as it relates to the elements g..{x). 

The part of the theorem that relates to the elements h..{^x) is proved in an 
entirely similar way. 

The solutions ^j,(x), • • • , g,^-{x) ^h^.(x), • • • , h^{x)'\ may be said to be 
the solutions associated with the determinant limits u^(^x), ^^.(a;) ■ • • [«;(«), 
v,j(x) ■■■'}. 

It is irwportayit to obtain explicit formulas for these functions g..(x), and 
we shall show that it is possible to obtain the solutions ^i(a;), ^3. (cc), • • •, 
g^j(x) for j =1, 2, •■•,n in terms of u.(^x), u..(x), •••, ?«i2...tci (fnd the 
operation X defined by the condition 

z <p(t)-j:<p{t)=4>(^)- 

t—x+l t—x 

The functions g'n ( cc ) , • • • > 6'„i ( ^ ) ^^re equal respectively to u^{x) , • - • , w^ ( x ) . 
It will accordingly be sufficient to obtain formulas for g^^ ( a^ ) » • • • ? g„^ ( cc ) in 



terms of the solutions </y(a;), • 
determinant limits. 

Form the determinant 


•' ffnj{^)U = '^^ 


2, 


• • • , ^ — 1 ) and 




ffni^^) 


9ui'") 




ffu(^) 


(23) d('=>{x) = 


^uC^' + l) 


9',2(^ + l) 




5'u(«' + l) 




g^^(x + k-l) 


g^^{x + k-l) 




g^^(x+k-l) 



If now in place of the elements gi-(x + 7?i ) ( »j > ) we insert their value as 
obtained from (1) 

n 

(24) g^.{x + m)= ^ a^^{x+ m-l)a^^{x + m-2) ■ ■ . a^g(x)ggj{x), 

T, <r,..., »=1 

in terms of the functions g^.[x), ■ ■ ■, g^.(x),v/e see that the above determinant 
may be expressed as a sum of terms each of which contains as factor a deter- 
minant which reduces to one of the determinant limits of order k . In fact we 
have obviously 

(25) e^"\x) = ^ ±_^[a,^{x)^±a,^{x+l)a,„(x)'^. . .]u^„„Xx). 

Hence this determinant is to be regarded as a known function. 

But the identity (23) may also be looked upon as a non-homogeneous linear 
difference equation of order k — 1 in g^^^x), and may be written 

(26) e('^-^Xx)g^Xx+k-l)+f(x)g^Xx+k-2)+ ■ ■ .±d<"-^Xx + l)g^^{x)=e<''\x). 
The associated reduced homogeneous equation obtained by equating the left 
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hand member to zero has the solutions [see (23)] gii{x), ■ • • , </i, ;,_, ( a; ), which 
are linearly independent unless 6'^''~^^(x) is zero. By a method analogous to that 
of variation of parameters for ordinary linear differential equations, let us put 



(27) 






and determine Iji^x'), l^(x), • • • , lJ^_^(x) by the equations 



(28) 



j:g,^{x+l)Al.{x) = 0, 
'g<7,.(a,+ 2)A^.(,^•) = 0, 



j = i 






6l«-i'(x) 

where as usual A(^ denotes the difference 4>(x -j- 1) — (j)(x). There are here 
k — 1 linear equations to determine the k — 1 unknowns Alj(x),- • • , A4_i {x), 
and their determinant is precisely ^*~''(x4- 1). Unless this determinant is 
identically zero, we have accordingly 

^l^j[X) - ^te-,)(a;)^«-l)(a, + 1)' 

in which W7^- j,_i(a!) is the cofactor of the element in the jth column and last row 
of the determinant ^*-"(a! + 1): 

(7u(a; + l) ■ ■ ■ ff^j_,(x + l) (j,j_^,{x+l) ■ ■ ■ ff,,„-i{x+l) 

ffn{^+^) ••• 6'i,i-i(« + 2) gij+,{x+2) •••</,,*_,(« + 2) 



(29) ± 



g^^{x + k-l) . . . ffjj_^{x+k-l)g,j^,{x+k-l) ■ ■ ■ cj,^,,_,{x+k-l) 
The resulting expression for g^^'x) ^^ 



(30) 






It remains to show that gyX^) ^^ ^^^^ defined is a solution of (26). We have 

7—1 

^i.(a') = Z5'i,(a')^,-(a')' 
5'u(«' + l) = £[(^u(«'+l)^,(«') + ^i,-(«'+l)A0(^)] 
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;,— 1 



g,,{x + k- 1)= £ \_{g,j{x + k -l)L{x) + g,j{x + k- l)A^(x)] 

= Z ,9,, {x + k-l)l.{x) + g^^^ 



Hence, if we multiply the last of these equations by 0*~'^(a;), the one which 
precedes it by y{r(x) and so on, and add corresponding members, the left hand 
member of the resulting equation will be precisely the left hand member of (26). 
Also all the terms except the last of the right hand member reduce to zero by 
reason of the fact that g„{x), • ■ ■, gi^i:-i{x) are solutions of the reduced equa- 
tion, and there remains only a term ^*^(a;). Therefore f/,,, (^) ^^ ^^^ desired 
solution : furthermore the formula (30) does not determine g'l,, (*) uniquely but 
only up to a sum of terms 

where J9, ( a; ) , • • • , p^_i (x) are arbitrary periodic multipliers since the operation 
^t=x ^® determined only up to an additive periodic term. But the result stated in 
Theorem II shows that g^j^ix) has precisely this degree of indeterminateness, and 
hence the most general determination oj" g^j.(^x) is furnished hy (30) unless 
&-^-'^\x) is identically zero, and this is an obvious case of failure. 

When a choice of g^X^) ^^^ been made, the remaining functions g^/^ix), ■ ■ ■, 
9nki^) which make up the ^th solution are determined. In fact the equations (24) 
for m = 1 , 2 , • • • , n form n linear equations in gr^ (a; ) , ■ ■ ■ , g'^(x) oi determinant 
not zero, since in the contrary case we should obtain a linear homogeneous relation 
between grj(x + 1), ■ ■ ■ , g^(^x + n), which contradicts the condition placed on 
gi{x) at the outset that it satisfies an equation of the nth order. By solving 
these equations we obtain 

(^2) g,{x) = (r.^{x)gXx + 1) + (r..Xx)g^{x + 2) + ... + <T.^{x)gXx + n) 

(i=l, 2, •••, n), 

where the functions o"; (a?) are rational in x. These equations determine 
S/zi^)^ ■ ■■,g^{x) in terms oi g^x). 

That none of functions ^^" ( a; ) , • • • , 0'"' ( x ) can be identically zero may be 
proved as follows : In the first place 0'*' (x) has a definite asymptotic form in 
the entire cc-plane, namely that obtained by replacing u^^__ (x) by s^^ (x) in 
the right hand member of (25). This is a consequence of Theorem I. But this 
form is 

^,{^) ••• %(•«) 

s,,{x+l) ... Sj,(x + 1) 

s,j(x + A; — 1) ••■ s^^{x+ k — 1) 
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as may be shown readily by means of (23). Recalling now the explicit form of 
the series S;.(x), we obtain 

(33) 6l(''>(a;)~a;''''^(/3je-'*)^. •• (/3,^e-*')^x'->+ •■• + '•'( <?,^ + -'^^ + ••• j. 

Since 0'^''\x) is represented asymptotically by this expression with respect to cc 
to the left of any line parallel to the axis of imaginaries, and with respect to v to 
the right of such a line (Theorem I), it is obvious that unless d,., d'-^, ■ ■ ■ are all 
zero, ^'"'(te) does not vanish at a great distance from the positive axis of reals. 
But these constants can not all vanish. Otherwise the above determinant 
vanishes (formally) and s^^^x), • • •■, f'uX^) ^^^ linearly dependent, and accord- 
ingly, by (32), s.^(x), • • •, s.,^(x) (i = 1, 2, . • •, w), are linearly dependent, in 
contradiction to (7). Thus (30) is a valid formula. 

§ 3. The Intermediate Associated Solutions. 

The application of the formulas (30) and (32) leads us immediately to a set 
of' intermediate associated solutions 

H'vA^)^ ■•■^9'„iX^) (fc = l,2, •••,n), 

analytic above and below lines A oo and B oo parallel to the axis of reals, and 
which have the property that 

with respect to x in any left halfj)lane and with respect to v in any right half 
plane. 
We take 

g[^(x) = ti,{x), ■.., g[^{x) = u^{x), 

and this first solution has the properties stated, by Theorem 1. We now define 
the solutions g[i{oa), ..., g'„oX^); g\A'^), ■•■, ff'„s(«')l •■■ in succession by 
choosing the summation operation S on (/)(x) to be defined as 

- cf>(x) - 4>{x + 1) 

and employing formulas (30) and (32). It will suffice to show that a solu- 
tion g' (x), • • • ) gl{x) for j = ^' will exist and have these properties under 
the hypothesis that the same is true for J =1,2, •••, k — 1. The formula 
(30), which serves to define 5''ij(x) if |jO,._i| > \p^\, is then 

^qzl^ ' ( \- V W .^ e''''>{x-]-v)m;^,_,(x + v) 

We have written m'^,^_-^{^x) to stand for what mj!^_i(x) becomes when g'.. 
replaces g.j . We shall first prove that the ^- — 1 series which appear here are 
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convergent. To this end we need to compute the asymptotic form of the typical 
term of the above series 

6"*-i'(a;)^"("-i'(a; + l)' 

and this is done by merely substituting the known asymptotic forms for the 
^'s as given in (33) and at the same time the known form of the elements 
mj ;j,_i . This latter [see (29)] is clearly given by 

(35) I „,(i) 1 

• • • (/3^_ie-")^a;'''+-+'V-i+'-^« +•■•+'■" j nij + — ■^ + . . . I 

as a result of our hypothesis. Combining (85) and (33) we find for the above 
term the asymptotic formula 

\ ) 0(*-i^(a;)0(*-'^(a;-|- 1) \pj {'^^ x^ 

The integer s. is zero unless (?j._j = 0. 

Since we have \pA > \p,\ for j < k except when | p^_^ | = | /3,J [see (8)] , the 
terms in (34) clearly diminish in approximately geometrical ratio for v large 
enough, and in consequence the series converge uniformly in the vicinity of such 
points. Thus g[,^{x) has been actually defined and is analytic above and below 
suitable lines Aoc and Boo . 

By the aid of relations (36) we may now determine the asymptotic form of 
each of the series in (34). In the first place, if each term of the jth series is 
replaced by its asymptotic formula, the series becomes 



Z ( - ) {x + vy^-'-'{x + vj' ( «. 

v=0 \Pj J \ 



^ X + V 



The error committed in the final bracketed series by breaking off at the ^th 
term is of order (.'*; + v)~'- if x lies to the left of the axis of imaginaries, and of 
order ii' if x lies to the right of this line, since the imaginary coefficient v of 
a; + I' is the same as the imaginary coefficient of x . These bracketed series are, 
of course, finite above and below suitable lines Aoo and Boo. 

Now separate the series into two parts : first the sum of the first m terms in 
which !< < J |a;| , secondly, the remaining terms. 

In this first set of terms we have 



{x + vy = xo(^i+^^+ ...^, 



where the error committed by breaking off the series at the ^th term is of the 
order (I'/a;)*. Hence the first m terms of the above series may be written 



262 G. D. BIRKHOFF : [April 

where Mj(^x) is limited. Now the terms involving M\x) are of the order a;~*, 



since 



!|:(^y-<^,(E)'Hfs(g)v 

As m becomes large we have also 



, M & fixed constant. 



%m 



T^ ' 



pji" 



pj - pk 

o all powers of x since the difference in absolute value is 



(£)" 



-vir 






< 



m' 



^ 



pj 



Hence the first m terms have the asymptotic form 
(37) (^Yx-^-^'x^'l 



(Xf 






with respect to x for x to the left of the axis of imaginaries, or similarly to the 
left of any parallel line. The terms after the mth diminish in approximately 
geometrical ratio for | v | sufficiently large, and aU contain a very small multiplier 
{pjjp.)'"'. Hence the second group of terms does not affect the above asymptotic 
form. 

Substituting the above expression in (34) and at the same time the known 
asymptotic formulas s^. (*•) for g[ . («), we find 



^;.(--«) ~ E ^•"■'^(P; e-" fx^' s, . + -f + 



i=^ 



o(1) 



or multiplying out, 



i$) 






(38) g\,Xx) ~ (P,e- fx^^x' { ^„ + f + • • 

with respect to x in any left half plane. 

By means of (32) we now obtain the asymptotic formulas for g'2^{x), • • ., 
g[,^{x). Thus we obtain n formal expansions which must of course be the con- 
stituent elements of some formal solution. The possibility that ^j,., I'j'^, • • • are 
all zero must be excluded since from that fact we should infer that all the coeffi- 
cients in the asymptotic expansion oi g'.^(^x)(^i = l,2, •••,w), and consequently 
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of ^*'(a;), were zero. The formal solution thus obtained must coincide with 
Sj^(x), ■ • ■ , s^^{x) at least if Pj, • • •■, p^ are distinct. For it is apparent that the 
only possible difference lies in a constant factor. By definition of ^*^(a;) in terms 
oi g^^{x) = g[^{x), • • • , g^^^ = g'^J^x), and in virtue of the fact that the asymp- 
totic form of 6^^ (x) has been shown to be that obtained by replacing these func- 
tions by Sjj(a;), • ■ ■, s^^(x), it follows that this factor must be 1. Hence we 
have in this case ^^j. = s^.,., ^!l = s'^y, • • •. Under our restriction, these equa- 
tions are algebraic identities persisting in the exceptional cases when p^ , ■ ■ ■, p^ 
are not distinct also. 

Thus we have demonstrated the existence of a ^th solution having the desired 
properties provided that \Pi^_i\ and \p^\ are not equal. 

We shall now indicate briefly the modifications that are necessary in the case 
I p^ _j I = I pj^ I . Suppose that p^ is one of a group of roots 

Pa) Pa+i^ ••••> P^ 
of equal absolute values. If p,, = p^, we have seen how to construct the ^th 
solution ; if p,^ = p^^^, one may invert the order of p^ and p^^, in (8) and thus 
construct a ( A^ + 1 )th solution possessing the requisite asymptotic properties, and 
so continue to the /8th solution. It is however necessary to know that the solu- 
tions thus constructed are associated with u^[x), u..(^x), ■ • ■ as stated. Accord- 
ing to the method used in § 2 there certainly exist ath, ( a + 1 )th, • • • , /9th solu- 
tions over congruent sets of points which satisfy all the limit determinant 
relations, obtained by interchanging the order of jO^, ■ ■ ■ , p^ in (8); for the alge- 
braic condition that this hold will be satisfied at all the approximations since 
these are actual determinants, and hence will be satisfied in the limit also. On 
this account at least certain determinations of the sum formula employed will 
give a ^th associated solution. Since however this ^th associated solution is 
determined up to a linear combination with periodic coefficients of the solutions 
before the ath, and the formula which gives it indicates a solution determined up 
to a linear combination with periodic coefficients of solutions before the ath, the 
general solution given by the sum formula is a ^th associated solution. Thus 
the solutions, obtained in this exceptional case as indicated, are associated with 
the determinant limits. 

It would be of considerable interest to determine the nature of these inter- 
mediate solutions near to the axis of reals. 

Considerations of symmetry show that there exists a set of solutions h' (x), 
• ■ • , h'^.(^x) (j = 1, 2, •••,«.) having the properties stated. 

§4. The First and Second Principal Matrix Solutions. 

We shall proceed now to prove the existence of certain two remarkable 
matrix solutions G(x) and H{x) of (4), and this will be done by means of the 
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formulas of § 2 and the intermediate solutions of § 3. It is clear that in order 
to define any matrix G(x) it is merely necessary to choose properly the opera- 
tions S(=j; which occur in (30). We shall use a contour integral to effect these 
summations, similar to that which Carmichael has used, writing 

<f){t)dt 



t=x J L 



,2n-AI'-l{a;-<) . 



e2^r-J{a;-0 — 1 ' 



X-1 




B' 



Fig. 2. 



where X is an integer and i is a path cc AB cc like that in fig. 2. The 
contour ^ i? is a simple curve which passes between x and x — 1 . The summa- 
tion function <f)(t) is supposed to be analytic 
except at points at a limited distance from 
the positive axis of reals and vyithin the 
region oo AB oc . Furthermore <^(<) is 
supposed to tend toward zero along lines 
A x , B CO sufficiently distant from the axis 
of reals, and in such manner as to make the 
above integral absolutely and uniformly 
convergent in the vicinity of any x not congruent to a point oi L . In 
this case the function represented by the integral will be analytic at a 
sufficient distance from the positive axis of reals, and wiU be independent 
of the particular lines A <x> , B oo or particular curve AB chosen. If we 
form the difference Ylt=x+\'i'{^) ~ 2(=:c^(0 ^*' ^^^^ ^^ given by an integral 
over ABB'A'A, the integrand being unchanged. But within this closed curve 
the integrand is analytic save for a pole at t = x; and by the calculus of residues 
this difference reduces to <^(^x). Hence this is a possible determination of 

Having made this preliminary study of the operation ^I^^^., we now proceed 
to the proof of 

Theorem III. There exists a matrix solution G{x)\_H (^x)'^ of (i) whose 
elements have the following properties : 
(a) the functions $';■(*) [^i(a')] satisfy the conditions 



9ii{^^) = Ui{x), 



^.■n(^) = '"i(^)' 



ffjii^) 9jA^) 



= u,j{x), 



= ^iji^)^ 



(6) the functions 9{{x)\^hi.{x)'\ are analytic throughout the finite plane except 
possibly for poles to the right (left) of and congruent to the poles of the 
elements of A(x)^A~^ (a;— 1)] ; 

(c) in any left (right) half plane g..(x)^h..(x)'j is asymptotically represented 
by s..(x) with respect to x. 
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Proof. The proof of this theorem will be made by actually obtaining the 
solutions yii(a;), •••, g^^ix); •••; g^S^), •••, 5i'„„(;c) in succession. The 
conditions of the theorem are evidently satisfied by taking the first solution to 
be 

g,,{x) = u,{x), ■ ■ ., g^,{x) = u^x). 

It suffices therefore to exhibit a solution gyX-'')i ' ' ' ■< 9 nki'^') which satisfies 
the conditions of the theorem under the hypothesis that the solutions g^.(x), 

• • ■, g^.(^x) {j <ik) have been constructed satisfying all of these conditions. 
This is the way in which the theorem will be demonstrated. 

The formula (30) gives the most general determination of gj^(x) which is 
possible if g^,X'x) is to be the first member of a ^th solution satisfying restric- 
tion ( a ) . We shall consider first the case | /3,,_j | > | p^ | . We choose the 
operation 5^,^^ as indicated above, and obtain 

C39') a (x\--Ta (x\ f .'i^y^i.i—K^-f) ^^''^ ( * ) "^A ^-^ ^ ^ ^ '^ . 

the integer X.j, will be chosen to be the least integer as great as the real pari of 
(40) ^ ■j^== ( log p^ - log pj ) . 

The functions g^jX^)-: ' ' '•> 9niX^) ^^® *^®^ obtained from (32). 

The first step is to establish the fact that the summation functions in (39) 
satisfy the restriction on the summation function 4>{t). Now 0'^''^(^x), 0'-''~^Xx) 
are analytic save at points to the right of and congruent to the poles of ele- 
ments a^.(^) ^y (25). Likewise, by the initial hypothesis and by definition of 
^j,ir—i(x) as the cof actor of the element in the Jth column and last row of 
^**~''(a; + l) [see the form (29)], the singularities of mjj,_i(x) (_/ = !, 2, 

• ■ -, k — 1 ) are restricted in the same way. Hence the summation function 

..,^ 0'''Ht)m,,,_,(t) 

fy^)- e('c-i)(^t)e("-'\t + i) 

will be analytic at all points sufficiently far from the positive axis of reals, pro- 
vided that 0'-''~^\t) does not vanish except at points within a certain distance of 
that axis. This fact is an immediate consequence of the relation (33) in which 
(fj., d'-^, ■ • • have been seen not all to be zero. It remains to show that the 
function 4>{t) tends to zero in such a way that the integrals appearing in (39) 
converge as stated for all positions of the lines A oo and i5 oo at a sufficient 
distance from the axis of the reals. In order to demonstrate this essential fact 
we make use of the intermediate solutions obtained in § 3. By Theorem II the 
general determination of the solutions g^.(x), • • •■> g„(^) associated with 
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the determinant limits is given by any one of them plus arbitrary periodic mul- 
tiples of the preceding ones. Hence we have 

(41) gij{^) = g'a{'->')]\M) + •■■ + 9'i,i-^{^)Pi-^j{^) + g'iA^) 

(i=.l,2, •••,H;i = l,2, -.fc-l), 

where the functions p^.(x) are periodic functions of period 1, analytic for x suffi- 
ciently distant from the axis of reals since the determinant |^!.(x)| ~ | s..(a;) | 
does not then vanish. Substitute these expressions (41), for i = 1 , in m^- ^_i(t). 
A reference to (29) makes it clear that the »Jj *_] (i) may be written as a sum 
of products of periodic multipliers analytic for x sufficiently distant from the 
axis of reals and of factors mj j._i(i), where nij,._^ signifies as before that 
which nij j._i becomes when g' is put in place of g . Hence the integrals will 
converge in the manner desired if the expressions 

tend to zero in exponential fashion along A oo and B oo . This is true if ^ oo 
and B oo are sufficiently distant from the axis of reals, in view of the asymptotic 
relationship (36) which holds with respect to v in the right half plane. 

The function g^^ix) is therefore analytic at a sufficient distance from the 
positive axis of reals, and by (32) the same may be said of the remaining ele- 
ments ^27f(^)' ■ ■ ■' 9nh{^) which constitute the kth solution. Working to the 
right across the plane by the aid of the original equations one proves readily 
that the only possible singularities of the functions of this solution will be poles 
congruent to the poles of the elements of A(^x) and to the right of them. 
Hence (&) is true for the 7otli solution. Also (a) holds in so far as it relates to 
this solution since every solution g^^x), • • •, g„^{x) determined by (30) had 
this property. 

Statement (c) has to do with the asymptotic nature of ^i,,(*)5 • • •» S'n/X**')- 
In order to prove (c) let us break up each integral in (89j into two others over 
contours L^ and X, as follows: the contour L^ is a fixed contour oo A^B^<xi 
(fig. 3). If X lies above A^qo, L^ consists of a loop-circuit to oo which includes 
the points x, x + 1, • • • within it but not x — 1, x — 2, • • • . If a; lies 
between A^ oc and B^oo, L.^ consists of a loop which includes x, x -\- 1, • • •, 
X + I within it, where x + ^ is the last of this series of points to the left of 
A^B^. If X lies below B^oo, L^ is a. loop-circuit of the same nature as when 
X lies above A^oo. Any combination of paths L^ and L^ is clearly equivalent 
to the single path oc AB oo of fig. 2. If x lies on ^ oo or i? oo or on a line 
congruent to a point A^B^ on the left, the tv^fo integrals are undefined. 

Let us consider first the part of the integrals contributed by L^ . The inte- 
grands are periodic functions of x of period 1 and therefore the integrals will 
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x-i-k 
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Fig. 3. 



also be periodic since the contour Z/^ is fixed. This part of the_;th one of these 
integrals may be written as either of the products 



o2n-A,iV/-la; 



^2n{\jk-i)V'- 



-Ix i g-2n-(A,v-l)V'-.( 






l(x-t) -^ 

dt 



It is apparent that in the first of these the second factor remains limited as v 
becomes positively infinite, inasmuch as e^" ^~'^ tends toward zero ; likewise that 
in the second of these the second factor remains limited as v becomes negatively 
infinite. Hence the part of ^j^ thus obtained is 

(42) Ci^=. gn(«')qyX^)+ ■■■ +g:,k-i{a^)qk-i.k{x), 

according as x lies above A^cc, between A^<x> and B^cc^ or below B^co. The 
functions 5'u(a;), •••, 5'*_i,*(a') are periodic and limited. This is obviously 
true unless x is in the vicinity of a point congruent to a point of the contour of 
integration. However it is true in this case also since each of the three expres- 
sions above is unchanged if the contour be varied continuously ; of course the 
function qf,.{x) is not one and the same analytic function in the three cases. 

Next let us consider the contribution (7j^ yielded by the integrals taken over 
Xj. Let us suppose first that x lies above A^ 00 or below B^ao . In this event 
the integrals in C^^ may be evaluated as a sum of residues, each of which may 
be expressed in two ways, as indicated in the following formulas 



(43) 
or 






y=i 



*-i " r Qik 



6l(")(x-|- v)m'j^^_^{x+ v) 



{x + v)6^''-^\x + 



\-_v)_-] 
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The first form is the obvious one derived from equation (39). If we observe 
however that we have, by definition of m^. j._i(.'c), 

gn{x + v+ !)• ••$']*_,(«;+ I' + 1) 

ffn{«^ + V + 2) . . . ff,^^_i{x + V + 2) 



k-l 



T,gij(x)mj^,,_,{x + v) = 



i=i 



9n{x) ■■■ffi.k-i{x) 

and then write in place of g^. their equals as given in (41), the resulting deter- 
minant may be simplified because the periodic coefficients are the same at 
x + v-\-l,x-i-v + 2, ••• as at x. By a combination of successive columns 
the determinant reduces to a new determinant identical with the above except 
that g' replaces g, i. e., reduces to 

'E.g[j{x)mj^^_,(x + v). 
i=i 

Thus the second of the above forms is obtained except that the two signs of 
summation will appear in inverted order. It is permissible to change this order 
if all the infinite series converge. This is certainly the case as the second form 
is precisely that which defined g'i,.(x) in § 3 [see (34)]. Thus we have for x 
above ^j oo or below B^^ oo 
(44) C^ = g[,(x). 

When X lies between A^ oo and B^ oo to the left of A^ B^ we use a form like 
(43). We will agree to denote (Tj, by g'u(x) for the sake of convenience, 
because then (44) holds in all cases and g[^{x) as thus defined will have the 
■asymptotic form s„. with respect to x in any left half plane, and with respect 
to V in any right half plane. In § 3 these facts have already been proved for 
g[^{x) as there defined, that is for x not between A^ oo and B^ oo . We need 
therefore to prove that g[^{x) ~ s^j^[x) with respect to x between ^^oo and 
B^oo only. The discussion is entirely similar to a discussion that has been given 
in § 3. In the first place the integrals in the expression g'i,Xx) for C^ may be 
■evaluated as a sum of residues 

where x+ I is the last point of the set ic, ic + 1 , • • • to the left oi A^B^. If 
this equation be compared with (34) one sees that the right hand members differ 
in that a finite sum replaces an infinite series. Proceeding precisely as we did 
with (34) we arrive at the asymptotic relation g'i^(x) ~ s^j^(x) with respect to 
X by means of the hypothesis made that gij(x), ■ • • , g„j(^x) are asymptotically 
represented by Sy(ic), ■ • ■ , s„j(x) with respect to x in any left half plane and 
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for _; = 1,2, • • • , k — 1. The sole modification consists in this that instead 
of having an infinite set of terms in which those after the mth (m the greatest 
integer less than J ] a; | ) do not affect the asymptotic form, one has only a finite 
set of terms in which however those after the mth do not affect the asymptotic 
form. 

The proof of (c) of the theorem in so far as it relates to ^'^(a;), • • •, 9„j,{x) 
may now be given. In the first place, for x above A^ oo, (42) and (44) give 

The dominating term will depend upon the dominating term of 

which begin the formal expansions. On dividing by the last of these, the 
exponents become 

The bracketed expressions here all have a positive real part less than 1 : this 
follows from the definition of Xj^ , • • • , \^_i ^ and the fact that | /)j._j | > | />,, | • 
Hence the coefficients of x are represented by points in the upper half plane 
below the line « = 27r. Therefore if x lies in the second quadrant, all of 
these exponents have a negative real part, and the last term dominates, that is 

9'u(a')~s,i(»') 

for X above A^ oo and to the left of the axis of imaginaries, or likewise of any 
parallel line. 

For X between A^ oo and £^ oo we have similarly 

9ik{^) ~ «n{^)9A^) + • • • + s,,i_i(a;)?*_i,i.(x) + Su:{x) 

and the last term clearly dominates again by (8). 
For X below S^oo, ff^^{'x) is represented by 

s„(a;)e^-(^»-')^=^-?,,(x) + • • • + s,,,_.(x)e^''(^---«^=^^^,_,,,(x) + s,,(x). 
The dominating term depends on the exponents 

.„.zi(,._, ,'^^^.)., . . .,2.vn(v...-i/^^^>, 0. 

The coefficients of x will in this case clearly lie below the real axis and above the 
line «) = 27r, and if x lies in the third quadrant, the last term s^j^(x) again 
dominates. 

Hence in any left half plane gi!,{x) is asymptotically represented by 5,j.(a;) 
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with respect to x, and of course g 21X^)1 '"'■> ffnki^) ^^^ ^* *^® same time 
represented by s^^^i^x) , ■ ■ ■ , s^^j^(x) respectively, as follows from (32). Thus a 
solution g^,,{x), • • •■< ff„j,{x) exists possessing properties (a),(6),(c) unless 
\Pk-i\ — \Pk\i ^^^^ ^^ unless /3j. is one of a group p^, p^^^, •••,/) of roots p of 
equal absolute value. The same argument which enabled us to treat this excep- 
tional case in § 3 may now be applied here ; that is, we may take p^ as the first 
of this group without violating (8) and construct a A;th solution which satisfies 
(a), (6), (c) corresponding to each of these values of p . 

The second set of solutions h^.{x), • • • , h^.(^x), (j = 1, 2, ■ ■ ■ , n) may be 
obtained by considerations of symmetry. 

There is one and only one solution (/f, J (a?), • • •■, ff^j^i^) \_^n{^)'< ''''^njC^)] 
(k = 1, 2, •••, n) analytic except for poles and asymptotically represented 
^y *!&(*)' ' ' ' ■< ^n7ci^) 6ven to terms of the first order in any left (right) half 
plane. The matrix G{x)\_Il(^x)'\ with elements S'i(»')[\ (if)] is termed 
the first (second) principal matrix solution of (4). 

In order to justify this definition it is necessary to show that there exists but 
one such solution g^^^x), •••, 5'„;,(^)[^i;j(^)' ''■' ^nt(^)] * Suppose if 
possible that there exists a second solution ^,j(cc), • • • , g„^{x) with the proper- 
ties stated. Now we can express gf,.{x) as follows : 

9a{^) = giii^)Pi{^) + •■■ +giS^)p.{^) (i=i, 2, -•.,»), 

where p^(x), •••■,p„{^) are periodic of period 1, inasmuch as ^i^(x), •••, 
^„y (x) (_;' = 1 , 2 , • • • , n) form a fundamental system of solutions. If we solve 
these equations for ^0^.(03), we shall have 

5'n(a') ••• 9ik{^) ■■■ gmi^) 



-(^)-i^ 



9„A^) 



9ui^) 



9nni^) 



9ln(^) 



9nl{^) ••■ 9nii^) ••• 9n,X^) 

Consider now some definite period-strip formed by two lines parallel to the 
axis of imaginaries and one unit apart, and let x become infinite in this 
strip. On account of the known asymptotic form of the elements g and g in 
the strip to terms of the first order we find at once, if i 4= ^) 

for X at either end of the strip. If we write z = e^" ''', this may be written 



p,{x) ~ s2T^/^l^'°^'«-'°^'"\log s)^*--*[0]. 
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The upper end of the strip corresponds to s = , and the lower end to z = oo . 
The above relation shows that p-{x) has a zero of at least the order X^,, at 
s = . Likewise, this relation shows that p^{x) has a pole of order \,, — 1 
at most at z— co. Consequently ^ ( a; ) s""''** is analytic at s = and vanishes 
at z= 00 . Moreover, by taking x far enough to the left, the original determi- 
nant expression for Pi{x) is analytic since every element of the two determinants 
is analytic, and the denominator is asymptotic to the formal determinant | s.. (a; ) | 
which is not zero if x lies far to the left. Hence p^( a;) »"''»*, for j =f= ^? is 
analytic for all finite values of z and vanishes at s = oo ; it must be identically 
zero. We conclude that 

Now let i be chosen so that s.^, the constant of the principal term of s.j^{x), is 
not zero, as is possible by (7). Consider any point Xg and the points x^ — 1 , 
CC|, — 2, . • . to its left. As x takes on these values, ffijX^)/ffu{^) tends toward 
1, and thus we see that p^,{x) is 1 at the point x^. That is, ^,j,(a;) and g'„X^) 
are one and the same function, in contradiction to our hypothesis. 

If the functions a..{x) had not been restricted to be rational, the only modi- 
fication in §§ 1-4 would be that in the statement and proof of (6), Theorem III, 
other singularities than poles would enter. 

§5. The Fundamental Periodic Functions. 

Either of the two principal matrix solutions G [x) and JI(x) of (4) may be 
expressed in terms of the other, e. g., 

(46) Q{x) = H{x)P{x), 

in which the elements Pij{x) of P{x) are periodic functions of period 1 which 
may be termed the fundamental periodic functions. We shall now determine 
the nature of these functions. The fact that will enable us to do this is that 
G(x) and J^{x) are composed of elements of known asymptotic form in any 
period strip which may be taken to be one and the same in both cases. 

On account of the gain in simplicity we shall make a preliminary transforma- 
tion. Let 

d(x) = (x — a)(x — ^) • • ■ (x — X) 

be the least common denominator of the set of rational functions a..(x), and 
write 

gi{x)=T{x-a)T{x-^)---T{x-\)gi{x) (<=1, 2, ■■•, »), 

where T(^x) is the gamma function. The transformed system will be altered in 
that the common denominator has been removed. The formal solutions of the 
new system will be those of the old multiplied by a certain factor, and hence 
will satisfy the restrictions imposed in § 1. 
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The new coefficients are polynomials in a;. It is then permissible to take all 
the functions a..(x) as polynomials, say of degree /tt or less, and this we shall do. 
The first principal matrix solution G{x) is constituted of elements analytic 
throughout the finite plane, in this case, inasmuch as the elements a^j ( a; ) have 
no poles in the finite plane. 

From the equation (46) we have 

(47) P{x)=^H-\x)G{x); 

but if we take our period strip far enough to the right, the elements of H{ x ) 
are analytic functions of x and the determinant of S{x) wiU not vanish in the 
strip inasmuch as its asymptotic form is given by the determinant | %(«)!• 
Consequently the elements Ptj(x) are functions of s = g^" '^ which are analytic 
except perhaps at 2 = and » = 00 , corresponding to the upper and lower end 
of the period strip. In the case when the functions a..{x)ia (1) are merely 
rational it is clear that the periodic functions Pij{x) are analytic or have poles 
«xcept at 2 = or s = 00 . 

It is necessary now to remove a certain ambiguity. The expression s..[x) 
represents the elements g'i(x) and h^\x) asymptotically in any right and left 
hand side of the plane respectively, and is a many-valued expression which is 
changed to 

(48) e^'"^^^^e2'"->^^s.^.(a;) 

after a complete positive circuit of a; = [see (6)]. Hence if the same deter- 
mination of s.j{x) is chosen for gij{x) and hij{x) along the direction of the 
positive axis of imaginaries, the elements ffij{x) and A^j(x) will have a different 
asymptotic form along the direction of the negative axis of imaginaries, since in 
one case we pass from the positive axis of imaginaries to the negative by a posi- 
tive partial circuit, and in the other by a negative partial circuit. (It must be 
remembered that g'^ ■ ( x ) and hfj(x) preserve their asymptotic form in any left 
and right half plane respectively.) Let the determination for h^j{x) be termed 
Sij(x). Then that of ffijix) is obviously obtained by a 'complete positive circuit 
and is given by (48) in the direction of negative imaginaries. 

In the upper half of the period strip gij{x) and Ih^i^x) have a common 
asymptotic representation s^j(^x), and this proves that 

P(a;)~ S-\x)S{x). 

But the elements of s»j(a5) and Sy(^x)oi S~^{x) and S{x) are of the respective 
asymptotic forms 

■^ 1 ( s§^ 
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[see (9) and (6)], so that in consequence the relation holds 
Therefore, if we think of ^^^(0;) as a function of s, we have 

The single-valued analytic function Pij{x) of z is accordingly infinite only to a 
finite order at » = and either has a pole there, or is analytic. For i = J it is 
«lear that Pij{x) is an entire function of z taking on the value 1 at » ^ , that is 

(49) p,,{x) = l+cll>e''"'-^''+ .... 
For i ^j we have obviously 

lim p^J ( 33 ) z~'"' = constant 

since X^ was defined as the least integer as great as the real part of 

log pj - log p, 
27ri/— 1 

Hence Pij{x) has a zero of order \j, at least, at s = , and may be written 

(50) p^j{x) = e^-^«^'=^-(4») + c^j)e2'^=^- +•••)' 

the coefficient of the exponential factor being an entire function of z. 

It remains to consider the functions Pij{x) in the lower half of the period 
strip. Here the functions gij{x) and A^.(a;) have the respective asymptotic 
forms 

e^'"'^-'''e''"''^-'Sij{x) and %(«), 

and the matrices G(x) and H(^x) are asymptotically represented by 

8{x) and S^{x) 
respectively where 

S^(x)= S{x){e'^>'^~"' e2'"-W=i"8.^. ) .* 

Hence we obtain 

P{x)-^ S-\x)S,{x) or \_S-\x)S{x)^{e^^'''^^^''e'-'-i-^^8ij) 

or pT" p^jX-'''+''i{e^'"''^-'^e^'"'i-'-^ ^ij ) . 

Thus each function Pij («) satisfies the condition 

limp,,(x). '-^^^r^^i -^== ) e--W-. = g,.. 



*The notation (S,;,) signifies a matrix with term Oij in the ith row andjth column. 
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We see again that Pij{x) becomes infinite only to a finite order at z = oo. 
This demonstrates th&t p..(x) is a polynomial in z and that Pij{x) (i =1=^') is 
equal to z'^ multiplied by a polynomial in s . 
In the case i =j we have 

so that the polynomial Pu{x) is of degree n'mz with leading coefficient e^'^''' ^' . 
In the case t 4=7 we have 

\\vcip^^(^x)z~'^~'^*^ = constant, 

2=aD 

SO that the polynomial coefficient of a** is of degree /^ — 1 . 

Thus we have proved the following theorem : 

Theorem IV. Ziet G(x) and ir(x) be the first and second principal 
matrix solutions of (4) in which the elements of H{ x ) are polynomials of 
degree /x at most in x , and write 

G{x) = H{x)P{x). 

The fundamental periodic functions p^^{^x) which constitute P{x) are in this 
case of the form 

p.^(x) = 1 + c^i^e'-^^-^" + • • • + e2'"-i-'-f e^-f'^^^ 

in which \j denotes the least integer as great as the real part of 

^;,r:^(logp,-log,,)- 

The simplest illustration is given by the difference equation for the Gamma 
function, when one has 

y(a;) = e'^^-'^-"^ !^^^~ , h{x)= T{x), p{x) = l - e'"''''^''. 
1 [I —X) 

A similar analysis of the nature of the functions pfj(x) at z= and 2 = 00 
may be made in the more general case when the functions a^j{x) are not 
restricted to be rational functions, and shows that p^■{x) is analytic or has a, 
pole at » = and z = 00 with leading terms of the same form as those of (51) at 
these points. If the a^ .(x) are rational, it is obvious that the functions Pij{x) are 
rational in e"" ^'^^ ; we shall not stop to investigate further the structure of 
these functions. 
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§6. The Principal Matrix Solutions in the Entire Plane. 

The results of the preceding paragraph make it possible to discuss the asymp- 
totic form of the elements of G{x) [and likewise of II{x)'\ in the entire com- 
plex plane. In the first place it has been shown that 

5'*{'^)~%(«') 

along any ray from a; = in the second or third quadrants. To investigate the 
asymptotic form along rays in the first and fourth quadrants we make use of the 
equation 

(52) G{x)=H{x)P{x) or g,j{x)=i,h,^{x)'p^^{x) (i,j = l, 2, ■■-,»), 

T=l 

in which the elements hfj(x) are of known asymptotic form %(a;) in any right 
half plane, and the elements Pij{x) are known functions of x . 

Consider now a ray from jc = in the first quadrant and lying above the axis 
of reals. Along this ray, by equations (51), we shall have 

p,^ (x)^ e'-Mj'i-^- cf, where X« = , c(«> = 1 . 

Thus the preceding equation for gij{x) leads to 

(53) g^{x) ~ s,, (a;)e^''M-M? + • • • + «,„(^)e^'''-^^^ci»\ 

a formula which determines the asymptotic form in the first quadrant. 

To determine the dominant term we need to compare the relative magni- 
tude of the terms 

„x:f,2iT\u^-ix f.X „x „2:TX„j-i -Ix 

If these be divided by pj, the logarithms become 

^'^^C^^^'^^«)^ »-••■• 

By the very definition of A,y , the coefficients e^p • • •, £^^ = 0, . • •, e„^- of x in these 
expressions will lie in the strip Q = v <i2ir. Let Pj , • . • , P^ be the points which 
represent these coefficients in the complex plane. On account of (8) these coef- 
ficients will have decreasing real parts and the points P^, • • • , P^ are ordered 
from right to left. As arg x increases from to Jtt the above logarithms main- 
tain their relative positions but are rotated about a; = through the same angle. 
The dominant term along the ray arg « ^ t is that one for which the corre- 
sponding point is farthest to the right (i. e., the logarithm has the greatest real 
part) when a rotation t has been effected. None of the terms after the j'-th 
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can dominate since the corresponding points are to the left of ^^ = for all 
values of T < 7r/2 . 

Draw the convex broken line P^P,, • ■ • P below which none of the points 
Pi, • • •, Pj\\Q and whose vertices are at some of these points. If the plane 
be rotated through an angle t, it is clear that these vertex-points, and only 
these, will correspond to dominant terms. Let the acute angles which the 
successive sides of P^ P„ ■ ■ ■ Pj make with the axis of reals be <l>^, <f>^, • • • . 
The critical rays along which the dominant term changes are then in order of 
increasing magnitude 

2~^i' 2"^"' ■■■■ 
For small values of arg x the dominant term is clearly given by 

s«(a;)e2'V=J-c(5', 
and as the first critical ray is passed this changes to 

Along the ray the form is given by a combination of two or more terms, a fact 
which we shall express by saying that the change is uniform along the ray. A 
similar change takes place along each of the successive critical rays until the last 
is passed when the asymptotic form becomes Sy{x). 

Likewise in the fourth quadrant the asymptotic nature is deduced from a 
relation 

g^^{x) ~ s.j(a;)e^''<''+^'^-»^^=^c',5-'^ + • • • + s<^(a;)e"^'"'^^»'e=''">'^-i" 

+ . . . + s,„ ( ,^ ) e 2.(^+x,.-i)^/-ix c(g-i) . 

The expressions s^^(x), • • • i s<„(a5) which appear here are the determinations 
for the functions h.^(x), ■ • -, h^^{x)'m the right half plane. Consequently the 
determination of s^ (x) obtained by a positive rotation from the first to the fourth 
quadrant along the negative axis of reals is this determination multiplied by the 
factor e^'^'^-'^e'''"''''"'. Hence in terms of the second determination 

ff^^{x) ~ s,,(aT)e''(^v-i)''=T-e-^'"-/^c$,^-« + • • • + s,j(x) 
After division the corresponding logarithms are 
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The coefficients of x are represented by points Q^, • • •, ^^. = 0, • • •, Q^^lir 
units below the corresponding points Pj , • • • , P^ (except for P. = ^^. = ). These 
points all lie in the strip — lit = ^j < and are given in order from right to left. 
As arg X decreases from to — 7r/2 , the figure formed by these points rotates 
through an angle — 7r/2, and the dominant terms correspond to the point farthest 
to the right. None of the terms after the_;th can dominate. 

In this ease we are led to the following rule: draw a convex broken line 
Q,\Q,i • • ■ 0,3 above which none of the points Qj , Q^ ' • • • ' 0,3 li® ^^^ whose 
vertices are at these points. Let the acute angle which the successive sides 
§1 §9 • • • Oi ™^ke with the axis of a; be ->/r, , ->/r^ , • • • . The critical rays 
along which the dominant term changes are then 

in order of decreasing magnitude. For small arg x the dominant term has the 
asymptotic form 

and as the first critical ray is passed, this changes uniformly to 

and similar changes take place along the remaining critical rays until the last 
ray is passed, when the form is Sy(a?). 

Theorem V. Let the coefficients aij{x) he polynomials of degree not 
greater than fi in x. Write 






where \j is the least integer as great as the real part of the negative of the 
first term in parenthesis, and where p^, p^, • ■ ■, p^ are the roots of the 
characteristic equation ordered according to decreasing absolute magnitude. 
Marh the two sets of points 

eij, e-ij, •••, ej_^^j, a7id e.^ — 1 , e.^j - 1 , •••, €_,._,, ^. — 1, 

and call them respectively 

P„ P„ ..., P, and §., Q„ ..., Qj. 

Construct the convex broken line P^P^ ••• Pj above which all the remain- 
ing points P lie, and likewise constt'uct the convex broken line Q, ^^ • • • Oj 
below which all the remaining ' points Q lie. Let the acute angles which 
the successive sides of P^ P ^ • • • Pj and §| Oe ' ' ' Oj ''^'^^^ with the 
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axis of reals be (j>^, (p^, • • • and -*]r^, -^jr^, • ■ • respectively. The critical rays 
between arg a,' = and arg x = 2'ir are then 

2" 9i' 2~'^'" ■■■' —2^^^' ~9, + yx^ 

in angular order. The asymptotic form of giji^x) changes uniformly from 

SiJx)e-''''oj''^''c<^i^] to Si^(x)e'''''fi/-^''cfj 
along the critical ray arg x = ^tt — (f>^_^ and likewise from 

Siy{x) e^-CA^ ;-i)''-i-'e-2'"-s ''-^c^,"-') to Si^{x) e-^'Cs/-''''^^^ e-^"-s''-^c<i^-'^ 

along the critical ray org x = — Jtt + i/f j . Between the last critical ray of 
the first set and the first of the last set the asymptotic form is given by s^^^x). 

It is interesting to observe that with g^.(^x^, '''■> ff „■(•*') t^'c^'C are associated 
at most j — 1 critical rays in the first and in the fourth quadrants so that, for 
example, when j = 1 there are no critical rays. But we have gi{x) = u^(^x), 
• • •, ^^(x) = «„(«) and this result is in accord with Theorem 1. 

There are two special cases of exceptional interest. The first is that in 
which p, , P2 ' ' ' '1 Pn ^'"^ °^ equal absolute value. The quantities e. . are in this 
case all pure imaginaries and the broken lines P^ , P^ , • • • and Q^, Qg, • • • 
reduce to segments of the imaginary axis, and we have (j)^ = 7r/2 , i/^j = 7r/2 , 
so that the critical rays fall along the positive axis of reals, i. e., do not exist 
since this line is excluded throughout. In this event the elements of G(x) 
preserve their asymptotic form throughout the entire plane, a result which 
agrees with that obtained in § 1 . The second case is that in which the argu- 
ments of Pj, ■ ■ -, p^ are equal. In this case the quantities e,y are all real. In 
the first quadrant the only critical ray is the positive axis of imaginaries. In 
the fourth quadrant one finds critical rays different from the negative axis of 
imaginaries. How shall we reconcile these unsymmetric and apparently contra- 
dictory results, since by Theorem III the positive axis of imaginaries is not a 
critical ray ? If one of the constants cW which determine the change in asymptotic 
form vanishes, the formulas are still valid but do not give the dominant term.* 
The true critical rays in the first quadrant are symmetrically placed with respect 
to those in the fourth quadrant. 

In general the effect of a zero value of cW (or of c'^"""'^) is to shift the 
actual change in asymptotic form to a secondary critical ray making a smaller 
angle with the positive axis of reals. The complete determination of its posi- 
tion may be effected when the constants c'-^ are known. 

The critical rays for these solutions have been determined by Galbrun by 

* Just as, for example, if u = x' we have « -^ e* ( ) along the positive axis of reals, but this 
does not give the dominant term. 
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use of the Laplace transformation. It may be noted that the above results hold 
when the functions a^j (x) are only restricted as in (2), since only the form of the 
functions Pij{x) at s = and z = oo has been made use of. An interesting dif- 
ference between the case when the functions a^ ( a; ) are rational and the case 
when they are not lies in the fact that there are a limited number of possible 
positions for the secondary rays in the first case, and an unlimited number in 
the second. 

§ 7. The Fundamental Pkoblem. 

The properties which we have thus far adduced completely characterise the 
solutions of a linear difference system (1) with rational coefficients, as is evi- 
denced by the following theorem : 

Theorem VI. Let G(x) and /I{x) be two matrices of single-valued func- 
tions, analytic save for poles in the finite plane and such that their elements 
have the property 

lim g^^ ( a; ) x-'^^py'' x"*"' = s,-,- and li m h^J(x) x~'^'' pj" x~''> = Sy , (| s^^. | =}= ) 

for X in any left or right half plane respectively, and furthermore let 

G{x) = H{x)P{x) 

where the elements of P{x) are fer iodic of period 1. Then G(x) and H[x) 
are the first and second principal matrix solutions of a system (4) in which 
the elements of Ai^x) are rational in x, at least if Pii P21 • • ••> P„ '^'^^ distinct. 
Proof. If we put 

A{x) = G{x+\)G-\x) = H{x+\)H-\x), 

as we may by the hypothesis of the theorem, it is clear that the elements afj(x) 
of ^ (a;) are single- valued and analytic in the finite plane except for poles. Fur- 
thermore, by direct computation, the elements ^^^- (a;) and A^(a;) of G~^{x) and 
II~^ (a;) are seen to be such that 

limg'^J(x)x>'''(p^e''yx'''= a-^j and lim Ay (a;)a;*''''(j0.e~*')*-"ic''' = cry, 

X=^ X=Xl 

in which a^ denotes as before the element in the tth row and ^'th column of 
the matrix inverse to (sy). The variable x lies in any left and right half 
plane respectively. Therefore we have 

n n 

(tij{x) = T.9ir{^ + 1)5't/(«) = J1K{x + 1)^'t;('«)' 

T=l T=l 

n 

= Z (x+lY^'+'\p^e--^Y+'ix + iyr{s,^ + €,^)x->^^{p^e-^)-'x-''r{a^j+ e\j), 

T=l 

Trans. Am. Math. Soc. 19 
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where e^^ and e^ tend to zero for all large values of x , that is, in overlapping 
right and left half planes. This may be written 

ay{x) = X>'e~'^ll+ ~] 'Ei^irPr'^rJ+e-rj), 

where e"^j tends to zero as x becomes infinite. 

Thus a^J{x) must have a pole of order /ti at a; = oo with leading coefficient 

n 

T = l 

and is a rational function of x . In order to find the roots of the characteristic 
equation (3) in this case, we make an initial linear transformation with constant 
coefficients on g^[x), • ■ ■, g^{x) in the linear difference system 

n 

£'i{^) = 'Eaij{x)gj{x) 
1=1 

which has been constructed. The constants s^ undergo the same transformation 
and may be reduced to S^j ; thus we shall have a^j = ^fjPj, so that the invariant 
characteristic equation is 

{P-Pl){p-P2) ••■{P-Pn) = ^ 

with distinct roots p^, • • ■, p„ not zero. Hence the above system satisfies the 
restrictions of § 1. According to the definition of principal matrix solutions, 
G{x) and H{x) must be the first and second principal matrix solutions, since 
the principal matrix solutions are the only ones which preserve a definite asymp- 
totic form s^j{x) even to terms of the first order in any left or right half plane 
respectively. 

It is evident that the w^ -f 2«, exponential constants Pt, r^, Sy together with 
the functions Pij{x) characterize the principal solutions more or less completely. 
Let us turn again to the case in which the coefficients afj(x) are polynomials of 
degree fi in x. In this case the equations (51) show that the functions Pij(x) 
involves n^/u, — n further polynomial constants c'^J. But among the exponential 
constants there are n which are not essential since we can alter any set Sj^., • • • , s^^, 
by an arbitrary multiplier and not affect the associated linear difference system 
(1). Hence there have been specified in all 

(^n^ -\- 2n) -}- (n^ fi — n) — n or w^(/ti-fl) 

characteristic constants counting only the ratios Sy^, s^^ , • • • , 5„j. as character- 
istic constants. This number is precisely the number of arbitrary coefficients 
in the n^ polynomials afj{x), each of degree /u, or less in x. Now the system (1) 
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determines the constants, but not uniquely since r^, • • • , r^ may be replaced by 
any determination r^ + 2m^ V—1, •••,»■„ + 2m^ i/— 1 ( wij , • • • , ot„ integers ) 
in the formal expansions. There is however exactly one set of characteristic 
constants for each such determination. If the characteristic constants are inde- 
pendent, one may conversely regard the coefficients in the polynomials a^^.( a; ) 
as functions of the characteristic constants. The following definition and theorem 
go to show that conversely the characteristic constants do describe all the essen- 
ial properties of the solutions : 

Definition. Two linear systems (1) are said to belong to the same class if 
one may be transformed into the other by a linear transformation with rational 
coefticients. 

Theorem VII. Two linear difference systems (1) in which Py, p^, • • • i ?„ 
are distinct and in which the functions a,y(x) are polynomials in x of degree 
fi will belong to the same class if their characteristic constants are the same. 

Proof. Let 

(54) 6y{x + l) = Ay{x)Gy{x) and G,{x + 1) = A^{x)G,{x) 

be two systems of this kind, expressed in matrix notation, which possess the 
same characteristic constants. For convenience we suppose both of the equa- 
tions transformed into a normal form in which 

(55) a,j{x) = pjiB.jX'^ + a^jjx--' +...)• 

This may always be effected by one and the same linear transformations with 
constant coefficients, 

G{x) = CG(x), 

namely by that transformation which reduces the constants Sfj to B^J , 

(S*,) = (e,,)(.,,). 

This is always possible by (7), and if one substitutes in the new formal expan- 
sions, it appears at once that the transformed elements of Aj^{x) and -i4.2(a;) 
have the form (55). The exponential constants pj , rj and the polynomial con- 
stants are of course unaltered by this transformation, and we have explicitly in 
each case 
(66) ^i = <^ (i = l, 2, •••,«). 

This simplification being effected, let us consider the relation between the first 
principal matrix solutions Gi(x) and G^(^x). Write 

(57) G,{x) = R{x)G^{x). 
We find then 

(58) R{x) = G,{x)G-'{x)=Hy{x)P{x)lH^{x)P{x)Y^=Hy{x)H-\x), 
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where H^{x) and II^(x) are the second principal matrix solutions and P{x) is 
the common matrix of fundamental periodic functions. It is now easy to dem- 
onstrate that Ri^x) is a matrix of rational functions of x. The elements aij{x) 
of this matrix have only poles in the finite plane and the relation (58) shows 
that 

(59) \im E{x) = S,j 

in the complete vicinity of x = oo on account of the common asymptotic form 
of (tj and G^, H^ and ^. In consequence of (59) the elements »'y(a;) must be 
rational functions of x, analytic at 35= oo and taking on the value h^^ at that point. 
Substitute the value of G.^(x) given in (67) in the first of the given difference 
systems. There results 

G^{x + \)=lR-^{x-V\)A^{x)Ii{x)-\G^{x). 

A comparison with the second of the given systems shows then that 

(60) A^{x)B{x) = R{x-[-\)A^{x). 

This is the condition that the two systems belong to the same class and are 
transformable, the first into the second, by (57). 

Even if this theorem is true there remains the possibility that not all of the char- 
acteristic constants are independent, although we are now certain that these 
constants effectively determine the transcendental properties of the solutions. 
We may however say something more of the dependence of solutions on the 
characteristic constants. 

In the first place we observe that the exponential constants are expressed in 
terms of the coefficients in (56) when the reduction to normal form has been 
effected. Now the principal matrix solutions of (1) were constructed by means 
of certain limit determinants and a process of summation based on a contour 
integral. All the processes employed are uniformly convergent in the coeffi- 
cients dfl and involve them analytically. It follows that the two principal solu- 
tions G{x) and If(x) of (1) are analytic functions of the constants a^j* as 
well as of X, and hence the elements of 

P{x)= G{x)H-\x) 

have the same property. This necessitates, of course, that the coefiicients c'^,'^ be 
analytic in the coefficients aV) . Let us take /a, , • ■ • , p^^, r^, ■ • ■ , r^ as fixed, 
which amoimts to choosing a*/? in (56) once for all. The coefficients c'/j are then 
single-valued analytic functions of the remaining coefficients a'-.'}, which appear 
in precisely equal numbers. 

Now only two things are possible: either the characteristic constants are inde- 
pendent and one may conversely regard the quantities a['J as analytic functions 
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of the quantities c'/j, or else they are dependent. In this latter case we should 
be able to choose 

af = ^nk ( ^ ) for all i,j,k, 

where (^;.j(X) is an analytic function of the parameter A., in such a way that the 
characteristic constants are independent of X. The functions a,..(x) are then 
dependent on \, so that we may write A[x)~ A(^x, X). Furthermore we can 
take generally 

A^{x) = A{x,\), A^{x) = A{x,\) 

as two determinations of A{x) which give equations (54) with the same charac- 
teristic constants, and we may take X ^ X^^ as a value of X for which the func- 
tions ^. .,^(ic) are analytic. The matrix R{x) is defined by (57), where G^{x) and 
G^{x) denote again the first principal solutions in (54). The elements '>'^^[x) 
are rational functions of x by the preceding theorem. 

We will now show that the functions r^^ (a;), which depend also on X, must be 
polynomials in x. To this end we make the further assumption that no two of the 
/It" roots of the determinant equation |^(ic, X")] = differ by integers. This 
is manifestly the general case. 

If one of the elements of li has a pole at a point, it will be possible to form 
a group of such points 

a, a+ 1, •• •, yS 

containing the special point and such that « — 1 and yS + 1 are poles of no ele- 
ment of M. Now if we consider ( 60 ) at a; = /3, we see that the right-hand 
member is analytic at this point, and that the left-hand member will have at 
least one element with a pole at a; = /3 unless l^j(a!)| = at x= /3. Hence 
we must have ^,(y8) = 0. Likewise if we consider the same equation at 
a; = a — 1 , we see that the left-hand member is analytic at this point, and that 
the right-hand member will have at least one element with a pole unless 
|^^(a;) I = at .T = oc — 1. Hence we must have ^2 ("' ~ ^ ) = ^ • But for 
X near to X^^, the roots of 

|^,(a;)|= |^(ic, XJI = and \A^{x)\ = A{x, \)\ = 

will nearly coincide and no two of them, as a — • 1 and /3 respectively, can be 
congruent and distinct on account of the hypothesis concerning the roots of 
\A(x,X^^)\ = . It follows that the elements of Ji{x) can have no poles in 
the finite plane, and hence are all polynomials in x. According to (59) these 
polynomials reduce to S^j at a; = 00 , and hence are identically equal to S,^ . 
The elements of A(^x, X) will not depend on the parameter X, contrary to 
hypothesis. 
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Hence the characteristic constants are independent, and the coefficients in the 
polynomials a^j^x) may be conversely regarded as analytic functions of these 
characteristic constants. A fundamental problem of the linear difference 
system is to study this functional relation and in particular to determine 
whether there exists a system (1) with coefficients polynomial in x of degree 
not greater than fi having any assigned set of characteristic constants. 
Princeton University. 



